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1 Introduction 


In 1950 W. Kuhn and A. Katchalsky produced a fiber which contracts in response 
to changes in pH (Kuhn 1950). In their demonstration a single fiber immersed in a 
fluid repeated raised and lowered a weight proportional to the hydrogen ion content 
of the solvent. It was proposed at the time that such a device could be used as a 
linear actuator. One of the major drawbacks however was the exceeding slow response 
time* on the order of minutes. Since then a number of innovations have made this 
artificial “muscle” worthy of further study. First, fibers have been produced which 
contract in seconds and even tenths of seconds (DeRossi 1987; Suzuki 1989). Second, 
some fibers support considerable loads, on the order of lOON/cm 2 (DeRossi 1988). 
The contraction rates and forces have become comparable, and in some cases even 
exceeding, that of a human muscle. Third, the physics of fiber contraction (or gels 
as they are more generally called) have become well understood (Tanaka 1973a - 
87b). Finally, technological innovations such as robotics and implantable artificial 
biological organs have created a demand for such devices. A number of contractile 
gel devices have already been constructed, such as a robot gripper (Caldwell 1990), a 
multifingered hand (Toyota 1990), and an artificial urethral sphincter (Chiarelli 1988; 
DeRossi 1985; 86). 

The design of these devices from a practical engineering perspective, particularly 
with regard to dynamic modeling and control has only been considered recently (Gen- 
uini 1990). 

In this paper we will consider the design of a simple linear actuator based on 
contractile gel fibers as well as the design and simulation of single mechanical linkage 
controlled by two antagonist muscles. The first section will discuss the design of 
a simple linear actuator and the second will present a dynamic model along with 
approximate parameter estimates. The third section will introduce a nonlinear sliding 
mode controller which achieves desired trajectory tracking for model inputs. 



Figure 1-1. Two antagonist artificial muscles control a single link. 

2 Design 

Although there are many possible actuator designs, this paper will focus on the 
direct chemical to mechanical energy conversion through l lie control of the hydro- 
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gen ion content in the muscle fibers. The basic concept is to use high performance 
miniature valves to the control the inflow of acid or base to modulate pH and control 
contraction. A schematic diagram of the proposed actuator is shown in figure 2-1. 
Two input lines containing 0.1M HGL and 0.1M NaOH enter the muscle at the base. 
A miniature value placed into the base material controls the inflow of fluid. An irri¬ 
gation system design to facilitate mixing is composed of numerous small tubes which 
are interspersed among the fibers. Finally, a drainage tube allows the waste fluid (i.e. 
salt water) to be removed. 

The tubes and wetting surfaces of the valves are constructed from teflon which 
is chemically inert and the fibers are made from Poly-Vinyl Alcohol PVA. Although 
many different contractile materials exist (Tanaka 1990), PVA has been used suc¬ 
cessfully by a number of researchers and has a relatively high tensile strength. The 
tendons are made of Spectra and the tendon to PVA connections are machined Del- 
ran. The fibers are affixed at their terminations with epoxy which is also chemically 
inert. 



Figure 2-1. Two fluid microvalues meter desired amounts of acid or base 
through an irrigation system to a bundle of contractile fibers. 

3 Model 

The dynamic model of the system is composed of three basic parts: the fluid 
conveyance system, the hydrogel contractile fibers and the mechanical linkage. The 
following sections will address each of these systems in detail. 

3.1 Fluid System 

A schematic diagram of the fluid system is shown in figure 3-1. ft is assumed the 
inlet fluid line is under a moderate pressure P 0 and a microvalve with controllable 
resistance R v modulates the inflow of fluid. The system of irrigation lines is modeled 
as a fluid resistance Rj and an inertance /. Finally, the compliant sheath into which 
the fluid flows is modeled as a capacitance C and the exit line, a resistance //,. 
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Figure 3-1. The fluid model consists of a pressure source P 0 at the inlet of 
the system. The valve, tubes and drain line are modeled as fluid resistances 
R v , Rt and R e respectively. Finally, a fluid inertance I is included to 
account for the mass of the fluid and a capacitance C to model the expansion 
of the actuator chamber. 


Valve 


Although there are commercially available metering valves which continuously modu¬ 
late flow, they are generally too large from this application. Commercial piezo-electric 
or solenoid fluid values also exist in the appropriate size, although they have only two 
states, open and closed. It may be possible to use binary state valves in a pulse width 
modulation scheme to meter fluid, but this may introduce undesirable affects such as 
water hammer, excessive part wear and slower response. Alternatively there are some 
experimental metering valves and pumps which are of the correct dimension (in fact 
some are significantly smaller than this application demands creating the possibility 
of very small linear actuators). It may also be possible to use molecular valves in 
the form of biological or artificial membranes whose porosity is controlled by small 
voltages. 

In any case it is assume some mechanical valve of the appropriate dimension can 
I ><'• constructed, either a commercial two-state valve or an experimental multi-state 
valve. Based on commercial data sheets for fluid resistance we can approximate the 
resistance of a valve by 

R v = a{bu v ) n , 

where R v is fluid resistance in N 1/,2 sm 4 , a = 4.09 x 10 3 N 1/,2 sm 4 , b = 2.5 x 10 4 m/v, 
and n = —1.7925 are constants. For two-state valves we will assume u v = w 0 + sgn( i<) 
and for multi-state valves u v = u 0 + u, where u 0 is some small value and u is the 
control signal input, 0 to 12 v. Thesi parameters where roughly based on values 
given for the LEE Interface Fluidic 2-way (normally closed) microvalve. This valve 
is energized with 12 v at 250mW, switches in 1.5 msec and operates in a 0 to 7 psi 
pressure rang®:' The relation of pressure to flow is given by 

R 2 Q 2 = p 

±t v Kq 6v 1 v-) 

where Q v is flow rate in m 3 /s, P v is pressure in N/m 2 and R v is fluid resistance 
N 1/,2 sm -4 . 
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Irrigation System 


Multiple porous tubes attempt to uniformly enter and mix the fluid to increase the 
response time of the actuator. Consider a single tube as shown in figure 3-2. The 
tube has an inner radius r = 0.76 mm and a straight length / = 10 cm after a 90° 
elbow of radius r e = 4r = 0.30 cm. Small holes of radius 0.1 mm are punctured in the 
teflon tubes. The holes are small enough so that an initial pressure P m is necessary 
to induce flow. This prevents back flow into the acid-base lines, but for this initial 
model we will assume P m is zero. The holes are place at 90° intervals about the cross 
section and every 1 mm along the length. 

Suppose the tube as we describe is non-porous and open at the end. If we assume 
a fully developed flow under an entrance pressure of P 0 =1 kPa (above atmospheric), 
the fluid flow is 


7T P 0 r 4 




8 pi ’ 


( 1 ) 


where n = pv = 10 -3 kg/m s (for water, 25°C), r = 7.60 X 10 -4 m, and / = 0.1 m. 
This yields a flow of Q = 1.3 X 10 -6 m 3 /s and an average fluid velocity of u avg = 
Q/irr 2 = 0.72 m/s. The Reynolds number is Re = vd/v = 1089, where v = u avg , 
d = 1.52 X10 -3 m and v = 10 _6 m _2 /s, which is below the critical value Rey c „- t = 2300; 
therefore we may assume the flow is laminar. For higher pressures, equation 1 yields 
Reynolds numbers in the turbulent range, but the significant resistance offered by the 
sealed porous tube reduce the flow rates and maintains laminar conditions. As we 
will later show the maximum flow rate is approximately Q m ax = 1-65 X 10 -7 m 3 /s, 
which yields u avg = 0.1 m/s and Re = 152, which is certainly laminar. Finally, the 
entrance region before the fully developed flow is approximately l e = 0.06 Re d = 1.4 
cm; therefore we can assume the boundary layer is fully develop within the entire 
length of the tube. 

The wall shear force will result in a head loss and thereby create non-uniformity 
flows through the small irrigation holes. However, since the flow will decrease along 
the length due to seepage, the head loss will not be linear as in the case of homogeneous 
circular tubes. To approximate the head loss we will assume continuous porosity and 
a uniform fully developed flow at each point x of tube. The wall shear stress is given 
by 


T w (x) 


4 pQ(x) 

7rr 3 ’ 


which results in a pressure gradient 


dP(x) —8 pQ(x) 
dx 7rr 4 
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If we assume a flow loss due to porosity 


dQ(x) 

dx 


—aP(x), 


we have the resulting differential equation 


dP(x ) 2 
d 2 x 


bP(x ), 


where b = Sap/wr 4 . If we solvfr this differential equation with boundary conditions 
P(0) = P 0 and Q{1) = 0, we have the head loss as a function of length 

where c = \/b = 3.4 m _1 . The head loss in the tube is not significant in this case as 
shown in figure 3-2. 



Figure 3-2. Head loss due to wall shear force in the circular porous tube 
can generally be ignored. 

For the purposef. of the dynamic model, it may be more more useful to repre¬ 
sent the entire porous tube system as a fluid resistance Rj. Following a simplified 
convention, we will assume the fluid resistance is given by 

R = Ro + pU 

where R is fluid resistance R 0 is the resistance due to the entrance, p is the resistance 
per unit length and l is the length of the tube. The entrance resistance is given by 

R 0 = ad n , (2) 
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where d is the inner diameter of the tube in meters, a = 4.1 X 10 3 and n = —1.8. For 
the tubes 0.152 cm in diameter R 0 = 4.6 X 10 8 N 1/,2 sm -4 . The resistance per unit 
length 

p = a p d np , 

where a p = 1.22 X 10 6 N 1/,2 sm -4 , n p = —2, and d is the inner diameter. For the tube 
described above, the resistance per unit length is p = 5.3 X 10 8 N 1/,2 sm -5 , thus the 
total resistance for the tube is 5.13 X 10 8 N 1/,2 sm -4 . Using equation ?? to compute 
the resistance of the irrigation holes, we have a resistance 4.0 X lO 10 N 1/,2 sm -4 for 
holes 0.1 mm in diameter. However with 80 holes/cm and 10 cm of tube, the total 
resistance is 1.6 X 10 8 N 1/,2 sm -4 . 

For the acid and base systems we will have six irrigation tubes each, resulting in 
a total fluid resistance of Rj = 9.65 X 10 7 N 1/,2 sm -4 The relation between pressure 
and flow is again given by, 

R\Q\ = P T . 


Inertance 


Although the inertance of the fluid is small, fast switching times may cause its effects 
to be important. The inertance is given by 


I 


p\_ 

A’ 


where p = 10 3 kg/m 2 is the fluid density, l = 0.1 m is the length of the tubes and 
A = 1.1 X 10 -5 m 2 is the total cross sectional area. For this system, the fluid inertance 
is I = 9.2 X 10 6 kg/m 4 and the dynamics are given by 


IQi = Pi- 

Suppose the fluid accelerates to 1.3 X 10 -6 m 3 /s in 10 ms under a pressure of 1 kPa. 
Then the inertance / would be 7.7 X 10 6 kg/m 4 . This is comparable to the value given 
for the model. So for switching times on the order of 1 ms, inertance effects become 
significant. 


Capacitance 


In this particular design, the fiber system is covered with a compliant sheath. Since it 
is flexible, it can be approximated as a fluid capacitance; that is, when fluid is forced 
into the chamber the sheath will stretch and resist the flow. If we model the sheath as 
a thin-walled uniform cylinder with closed ends, the elemental stress tensor is given 
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by o r = 0, ere = Pr/t, cr z = Pr/2t and = 0, where P is the pressure, t is the wall 
thickness and r is the radius. Given the generalized stress-strain relations, 


e 


r 


£6 




^{(Te + (T Z ) 
-^((70 - va z ) 

^(a z - vag), 


and the strain-displacement equations 


A r = reg 
A / = le z , 


we can solve for the change in volume as a function of pressure 


AV 


7rr 3 /(3 — 4z/)\ 
2Et j ' 


where E is the Elastic Modulus and v is Poisson’s ratio. Therefore the fluid capaci¬ 
tance is 

n _ 7rr 3 /(3 — 4z/) 

“ 2Et ' 

If r = 1 cm, / = 10 cm, E = 10 6 N/m 2 (approximate for rubber), v = 0.5 and t = 0.1 
mm, 

C = 1.6 x 10“ 15 m 5 /N. 

With a fluid flow of Q = 1.3 X 10 -6 m 3 /s, this capacitance will become significant 
if the pressure changes more than 1 MPa in a 1 ms. Therefore we will ignore the 
capacitance. 

Drain 

The drain line consists of a single tube d e = 0.3 cm and a check valve. Employing 
equation ??, this results in a fluid resistance of R e = 2.22 X 10 8 N 1/,2 sm -4 , for a drain 
line 20 cm in length. 

Fluid dynamic equations 
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The elemental pressures are related by 




and the flows by 


Q — Qv — Qi — Qt — Q e* 


Substituting the elemental relations, we have 


Q = \ [Po - (fit + 4 + H‘)Q 2 } 


where Q is the flow into the fibers. Note the steady state flow rate for the open valve 
is _ 


Q ss 


Po 

Rl + R 2 T + Rl' 


For a fully open value Q ss 
is 


1.65 X 10 7 m 3 /s. Note the total volume of injected fluid 


V = Q. 


Mixing and Diffusion 

When the fluid (either acid or base) is injected into the fiber network there is a certain 
lag time before the entire system uniformly equilibrates to the new concentration. We 
could model this mixing as the diffusion of HC1 or NaOH in water. If the tubes are 
uniformly spaced as in figure 3-3, we can draw concentric circles around each with 
radius r 0 equal to half the distance between the tubes. The mixing time could be 
approximated, as the time for the fibers at radius r 0 to experience a certain concen¬ 
tration of fluid. As a simplification we can approximate the diffusion with the one 
dimension diffusion equations 


J 

dc 

dt 


d 2 c 
dr 2 ’ 


where J is the diffusion flux per unit area, c is the concentration in moles per volume 
and D is a constant. The solution is the standard diffusion model 


c = c 0 1 
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where 


erf(x) = 



e- y 2 dy. 


The boundary concentration c 0 , of course, decreases as particles diffuse into the liquid. 
However, as a first order approximation to the diffusion time, we assume c 0 is constant. 
Then the concentration at r 0 is one-half its initial value when t = 1.4 111 s for D = 
3.05 x 10 -4 m 2 /s (0.1M HGL a§ 25°G) and r 0 = 1.5 x 10 -3 m. 



Figure 3-3. Arrangement of acid and base lines in the cross section of the 
muscle. 

The diffusion model, however, is highly inaccurate since the exit velocity from the 
pores in the teflon tubes is relatively high and together with the motion of the fibers 
will cause turbulent mixing, not diffusion. Therefore as a simple approximation, we 
will assume the fluid mixing and diffusion is a first order lag, hence 

T m c + c = c 0 , (3 ) 

where c is the concentration of acid or base in the solvent, r m is a time constant 
which we will take to be r m = 1.4 111 s (based roughly on the previous discussion) and 
c 0 = Vi/Vj ratio of solute to solvent, where Vr is the total volume of the actuator, 
V T = 2.92 x 1(T 5 m 3 . 

3.2 Contractile Fibers 

pH 

The contractile fibers respond roughly linearly to pH, so it is important to calculate 
the pH in the solution based on the amount of acid or base entered into the chamber. 
An acid is a chemical species having the tendency to lose a proton and a base is a 
species tending to accept a proton. The quantitative measure of acid strength is the 
acid dissociation constant; that is the equilibrium constant for the reaction 

HA + H 2 0 = H+O + A- A' = 
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Along with this reaction there is always the self dissociation of water 


[K+0][0H~] = K w . 

If the original amount of acid A added to the solution is [HA] 0 , then the material 
balance relation is given by 


[HA] 0 = [HA]+[A"]. 

Also the solution must be electrically neutral, hence 

[H+O] = [A-] + [OH-]. 

Solving these equations for the hydrogen ion concentration yields the cubic 

[H+0] ([H+0] - A-„/[H+0]) 

[HA]„ - [HJO] + A'„/[HjO] 

However, we may assume that the Hg 0 from water is negligible compared to that 
from the acid, hence [Hg 0] = [A - ]. We may assume that for strong acids there is 
complete dissociation; that is [A - ] = [HA] 0 , which yields the trivial result 

[HA], = [H+0]. 


Similarly, for base a B we have equilibrium 


B + H 2 0 =+ HB +“ OH K 


[HB+]["0H] 


mass balance, 


[B]o = [B] + [ + BH] 


and electric neutrality 

[BH+] + [H+0] = [“OH]. 

Now we can assume total dissociation 


[B] 0 = [ + BH], 

which is a valid except for very high concentrations. Also we can assume the contri¬ 
bution of “OH from self-ionization of water is negligible, hence 

K K 

3 J = FMJ = lBp 
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Considering the strong acid and base together we can assume 


[H+0] 


[HA]„-[B]„ [HA]„ - [B]„ > 0 

K W /([B\ 0 — [HA]) otherwise 


Since pH is pH = — logfH+O] , we have 


-log([HA] 0 -[B] 0 ) [HA] 0 — [B] 0 > 0 

14 + log([B] 0 — [HA]) otherwise 


Diffusion 


Besides the lag due to mixing among the fibers, there is a lag due to the diffusion 
within the fibers themselves. This lag is due to bulk macromolecular rearrangement 
of the polymers in the network and is approximately proportional to the square of 
the cross sectional diameter. The time constant is 


r = cd 2 , 

where d is the diameter in meters and c = 2 X 10 9 s/m 2 . Thus for fibers 10 /mr in 
diameter contract in approximately 0.1 s. This can also be approximated to the first 
order by 

rpH + pH = pH s , 

where r = 0.2, pH is the pH in the fiber and pH s is that in the solvent. Notice this 
time constant is substantially slower than the one given in equation ??. Therefore 
we will ignore the lag due to mixing and consider only the diffusion within the fiber 
as the major delay. 

Hydrogel Physics 


A gel is a cross-linked network of polymers immersed in a fluid, which can undergo 
a volume phase transition in response to changes in external conditions such as tem¬ 
perature, solvent, pH, light and electric fields. Three competing forces act on the 
polymer gel network: rubber elasticity, polymer-polymer affinity and hydrogen ion 
pressure. Competition between these forces, collectively called the osmotic pressure , 
determine the equilibrium volume of the network. The phase transition has been an¬ 
alyzed in terms of the mean field theory of swelling equilibrium of gels. The osmotic 
pressure of a gel is given by Flory’s equation, 


7T = —- 


NkT 

r, , . ,, A F 


■ <f> u^x 

Vl 

+ ln(l ^)+2kT^\ 

+ vkT 

.2 <t>0 /. 


+ vfkT -f 
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where N is Avogadro’s number, k is the Boltzmann constant, A F = AH — TAS 
represents the difference between the free energies of a polymer segment-segment and 
polymer-solvent interaction, and AH and AS are the enthalpy and entropy respec¬ 
tively. In equilibrium, the osmotic pressure is zero. Thus by solving the above equa¬ 
tion, we can theoretically compute the fiber force. However because of its complexity 
and because of the availability of empirical data, we can describe the contraction of 
the fiber with the simple model shown in figure 3-4. The intrinsic force per unit area 
is f m is approximately 


fm = 1-0 x 10 5 N/m 2 (7 — pHj), 

assuming a total cross sectional area 2.92 X 10 -4 m 2 , the intrinsic force is F m = 29.2 
N(7-pHj). Also the intrinsic stiffness is 

k m = 1.43 x 10 3 N/m. 

and viscous damping is roughly estimated as 

b m = 30Ns/m. 





Figure 3-4. Simple model of contractile fiber. 

3.3 Mechanical System 

A model of the mechanical system is shown in figure 3-5. Assuming linear com¬ 
ponents for this part of the model, we have 

JO + (B + 2 R 2 b m )0 + 2 H 2 k m 9 = R(F 1 - F 2 ), 

where 
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J 

B 

R 

bm 

h 

n m 

e 

Fi 


1.28 x 10 -4 kg m 2 
neg. 

1.7 cm 
0.7 N s/m. 

1.43 x 10 3 N/m 
rads 

29.2 N(7 — pH•) 


Link inertia 
Joint damping 
Pulley radius 
Actuator damping 
Actuator stiffness 
Joint angle 

Controllable force for i th actuator 


I 


I 




Figure 3-5. Single link controlled by two antagonist muscles. 

3.4 Summary 

Using the model as outlined above we have the following system ecpiations 

Qi = ( P 0 -(R] + Rl)Q]-a\b{u 0 +\u\)) 2n Q])II 
Vi = QiSgn{u) 

P Hi = (-pH.-log(10- 7 + U/U r )/r 

6 = (-(B + 2R 2 b m )6-2R 2 k m 6 + R(F 1 - F 2 ))/J. 


However before we design a controller for this system let us consider a number of 
simplifications. First, the lag due to the inertance of the fluid is negligible compared 
with the dynamics of the rest of the system. This was verified in a full dynamic 
simulation of the system, although it is fairly obvious since the time constant of the 
fluid flow is on the order of 10 ms while the diffusion lag is about 100 ms. This 
assumption allows us to compute the flow Q directly from the input u. Thus we 
can reduce the order of the system by two and take the input to be the flow rate 
Qi. In order simplify this problem further let us assume a single-input single-output 
system, by taking the co-contraction of the actuators as an offset force and controlling 
position through feedback through a single actuator. With these initial assumptions 
let us write the system ecpiations in the form 

X = f(x, u). 
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Specifically, 


x\ = u 

x 2 = ai x 2 + a 4 log(a 2 + a 3 X 1 ) 
x 3 = X 4 

X 4 = a 4 x 4 + a 5 x 3 + a 6 + a 7 X 2 

where 

x\ = V(cm 3 ) 
x 2 = pH 
x 3 = 6 

X4 = 6 


and 


a 4 = -33 s 1 

a 2 = 10“ 7 

a 3 = 0.343 cm -3 

a 4 = -135 s _1 

a 5 = -6457 s -2 
a 6 = 2.89 x 10 4 s' 2 

a 7 = -3.87 x 10 3 s“ 2 


4 Control 

We wish to control the joint angle x 3 = 6 using the simplified model presented 
above. We would also like the controller to maintain zero tracking error in the presence 
of modeling errors. Therefore, we will design a non-linear sliding mode controller. 

However before we attempt a construct a controller for this fourth order system, 
let us initial consider one further simplification. Based on the numerical estimates of 
the model parameters the damping ratio of the linear mechanical dynamics is ( = 1.4. 
Therefore let us initially assume the inertia can be neglected, and produce a controller 
based on the reduce third order system, 

x\ = u 

x 2 = ai x 2 + a 4 log(a 2 + 03 ^ 1 ) 
x 3 = a 4 x 3 + a 5 + a 6 x 2 
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where 


and 


x 1 = V(cm 3 ) 
x 2 = pH 
x 3 = 6 


a\ = -33 s _1 

a 2 = 10 -7 

a 3 = 0.343 cm -3 

a 4 = -48 
a 5 = 201 

a 6 = -28.6 


Now the controllability canonical form , 

x ( n ) — f(xf + &(x)rt, 

where x = [ x, x,. . . , , can be derived by differentiating the equation for the 

joint angle three times. 


x 3 = a 4 a 6 log(a 2 + a 3 X\) + (a 4 a 6 + a 4 a 6 )x 2 + ofx 3 + a 5 a 4 
= a 4 a 6 (ai + a 4 )log(a 2 + a 3 x 4 ) + (a 4 a 6 (ai + a 4 ) + a\a G )x 2 

CL\(Iq 

+ 03^1 

Let x = x — Xd be the tracking error on x and define the sliding surface to be 
S(t) = {x £ 3? 5 |s(x,t) = 0}, where 

d 
dt 

for some positive A. Hence 

s = x + 2Xx + X 2 x. 

Thus maintaining zero tracking error is equivalent to the first-order stabilization prob¬ 
lem in s. In addition a bound on s implies a bound on the tracking error vector; that 
is if \s(t)\ < $, then |5i(t)W| < (2A)*e. We must now design a controller such that for 
all x ^ S(t), 


s(x,t) 



^ 4^3 + Cl4 a 5 + 


(a 2 
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for some positive rj. The estimation of / is given by / and the bounds on the error is 

I/-/I < F. 


Similarly, the term b is estimated by b — V^min^max where 0 < b ^ n < b < 6 max . If we 
let f3 — xj ^max/hiin? then. 


and the control law 
for s given above, 


u = b 1 (u — ksgn(sj) 
k = f3(F + v) + (f3 - l)|u|, 


and 

u = —/ + — Xx — X 2 x. 

In this case we chose A = 10 s -1 to be approximately equal to the natural frequency 
of the system while the estimate errors to be approximately 10% the nominal values. 


5 Simulation 

In order to get a sense of the model dynamics a 1Hz sinusoidal input was given to 
the system. A plot of the pH as a function of time is shown in figure 5-1 and the joint 
angle in figure 5-2. The sliding mode controller discussed in the previous section was 
implemented and simulated for a desired joint trajectory 

Xd = sin(27rt) + 1.5. 

Figure 5-3 shows the pH and figure 5-4 shows the desired and actual joint position. 
Notice desired tracking is achieved after approximately 0.5 seconds. Figure 5-5 shows 
the effect of neglected link inertia. Note the joint trajectory is offset by approximately 
10° . Fourth and fifth order sliding mode controllers were attempted to compensate 
for the link inertia and fluid inertance, but were found to be numerically unstable 
in the calculation of the higher order derivatives. Numerical filters corrected this 
problem, though the simulation times became unacceptable long. A more judicious 
choice of numerical parameters could possibly correct this problem. 
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Figure 5-1. pH in the polymers fibers in response to a 1Hz input signal. 
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Figure 5-2. Joint angle. 
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Figure 5-3. Response of pH to control of joint trajectory. 
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Figure 5-4. Desired and actual (modeled) joint trajectories. Zero tracking 
error is achieved after one second. 
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Figure 5-5. Effect of neglected joint inertia. 


6 Conclusion 

A sliding mode controller produced reasonable tracking performance for input sine 
waves on the order of 1-2 Hz, but deviated significantly at higher frequencies. This 
seems reasonable based on the relatively long lag times introduced by the polymer 
diffusion. It has been possible more recently, to produce even smaller fibers (about 1 
/mi), while using UV induced transectional cross-linking to increase strength. These 
and other innovations should allow faster response times for practical actuator devel¬ 
opment and controller design. 
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